N 72- m?7 


NASA CONTRACTOR 
REPORT 


CO 

CO 

CM 

CM 



<c 

on 


<C 


CASE FILE 

COPY 


DRAG COEFFICIENTS FOR SPHERES 
IN FREE MOLECULAR FLOW IN O 
AT SATELLITE VELOCITIES 


by John W. Boring and Robert R. Humphris 

Prepared by 

UNIVERSITY OF VIRGINIA 
Charlottesville, Va. 2 2901 
for Langley Research Center 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • MARCH 1973 



2. Government Accession No. 


1. Report No. 

NASA CR-2233 


4. Title and Subtitle 


DRAG COEFFICIENTS FOR SPHERIS IN FREE MOLECULAR FLOW IN 0 
AT SATELLITE VELOCITIES 


7. Author(s) 

John W. Boring and Robert R. Humphris 


9. Performing Organization Name and Address 
University of Virginia 

Research Laboratories for the Engineering Sciences 
Charlottesville, VA 22901 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 205*r6 


3. Recipient's Catalog No. 


5. Report Date 
March 1973 


6. Performing Organization Code 


8. Performing Organization Report No. 


10. Work Unit No. 
863 - 11 - 00-01 


11. Contract or Grant No. 
HASl-2538 


13. Type of Report and Period Covered 
Contractor Report 


14. Sponsoring Agency Code 



16. Abstract 

The drag coefficients for the Echo I and Explorer XXIV spherical surfaces in sm 0 
environment have been experimentally determined over an energy range of U to 200 eV. The 
experiment was performed by generating a beam of atomic oxygen ions of the proper energy, 
neutralizing a portion of the beam, and then allowing only the neutral 0 particles to 
strike a very sensitive torsion balance. The momentum transferred to the surface was 
determined from the deflection of the torsion balance. At the lover energies, the more 
intense ion beam had to be used instead of the neutral beam. The drag coefficients are 
found to be slightly greater than 2 at energies corresponding to satellite velocities. 


1 7. Key Words (Suggested by Author(s) ) 
Drag coefficients 
Atomic oxygen 
Echo I and Explorer 2k 
Momentum transfer 


19. Security Classif. (of this report! 
Unclassified 



For sale by the National Technical Information Service, Springfield, Virginia 22151 


22. Price* 

$ 3.00 

























TABLE OF CONTENTS 


Page 

List of Illustrations v 

SECTION I INTRODUCTION 1 

SECTION II EXPERIMENTAL METHOD 4 

A. Equipment 4 

B. D.C. Method (Neutrals) 14 

C. Swinging Method (Neutrals) 16 

D. D.C. Method (Ions) 17 

E. Surface Condition 20 

SECTION III RESULTS 22 

SECTION IV CONCLUSIONS 34 


REFERENCES 


35 



LIST OF ILLUSTRATIONS 


Page 

Figure 1 Experimental System 5 

Figure 2 Typical Mass Spectrograph 7 

Figure 3 Torsion Balance 9 

Figure 4 Typical Data Traces 18 

Figure 5 Photomicrograph of Surfaces 2 3 

Figure 6 Momentum Transfer Results for Echo I ... 24 

Figure 7 Momentum Transfer Results for 

Explorer XXIV (Lines Vertical) 25 

Figure 8 Momentum Transfer Results for 

Explorer XXIV (Lines Horizontal) 26 

Figure 9 Typical Results as a Function of cos9. . . 28 

Figure 10 Drag Coefficients for Echo I 30 

. Figure 11 Drag Coefficient for Explorer XXIV .... 31 

Figure 12 Average Drag Coefficients 32 

Figure 13 Side View of Smooth Surface and Surface 

with Grooves Horizontal 3 3 


v 



DRAG COEFFICIENTS FOR SPHERES IN FREE 
MOLECULAR FLOW IN 0 AT SATELLITE VELOCITIES 

BY JOHN W. BORING AND ROBERT R. HUMPHRIS 
UNIVERSITY OF VIRGINIA 

SECTION I 
INTRODUCTION 

One method of investigating the density of the earth's upper 
atmosphere is by observing the orbital decay of satellites. From 
a knowledge of the orbital parameters as a function of time one 
can infer the drag force produced by the atmosphere and hence it 
is possible to calculate the atmosphere density, if the drag 
coefficient is known. For the free-molecular flow conditions 
that exist for satellites in the upper atmosphere, a knowledge of 
the drag coefficient is equivalent to knowing the atmospheric com 
position and the manner in which individual molecules exchange 
momentum with the satellite surface. There is, however, very 
little experimental information available concerning molecule- 
surface interactions for relative velocities in the satellite 
range, and a consequent uncertainty in any corresponding theoreti 
cal calculation. 

In April 1963 the University of Virginia began a program 
to study in the laboratory the transfer of momentum from 
atmospheric molecules to solid surfaces, especially for 
relative velocities in the satellite range and for surfaces 
that are samples of actual satellite material. The program 
has been sponsored by NASA Langley Research Center under Contract 
NASl-2538. The original proposal [1] described the general ap- 
proach used, and earlier reports have described the system de- 
veloped to make the N 2 measurements [2] , preliminary results for 
N 2 molecules on several surfaces [3] , a study of various methods 
for producing a monatomic oxygen beam [4] , and the final 
momentum transfer results for N 2 on surfaces of Echo I, 

Echo II, and unpainted Explorer XIX satellite material [5]. 

The present report describes the final momentum transfer 
results for atomic oxygen, 0, on surfaces of Echo I and 
unpainted Explorer XXIV satellite material. Since much of 
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the experimental apparatus and procedure is similar to that 
presented in the final ^ momentum transfer results report [5] 
some of the sections of that report are repeated here. 

It was originally intended that measurements would also 
be made for several satellite surfaces which were painted with 
epoxy paint, but it was found that these electrically insulat- 
ing surfaces could not successfully be employed with the 
present technique, since small electrical charges on these 
surfaces produced forces that tended to obscure the small 
forces that were to be measured. A brief discussion of 
estimated drag coefficients for the painted surfaces is given 
in the final section of this report. 

Consider now the drag on a satellite moving through a 
rarefied atmosphere (free molecular flow) with a speed large 
compared to the thermal motion of the atmospheric molecules. 
The drag coefficient can then be expressed as [3] 

p 

C D = J — 2 = 2tl + i / s cos6da] (1) 

2 Apv o ° 

where F is the drag force, A is the cross-sectional area of 
the satellite projected on the plane normal to the direction 
of motion, p is the atmospheric density, v q is the relative 
velocity of the satellite through the atmosphere, P Q is the 
corresponding molecular momentum, P m is the average component 
of momentum of reflected molecules along the direction of 
motion (taken positive when opposite to P q ), 9 is the angle of 
incidence of molecules (measured from the normal to the 
surface) striking an element of surface da, and the integral 
extends over the surface of the satellite. The above expres- 
sion may therefore be used to calculate the drag coefficient 
-fo r — a — bo - 4 - y — o - f — e on-ve-x— sh-ape — (-s o — t - h at — dou-bie— pe^lecti-ono — are — not. 
possible) moving through a one-component atmosphere if. the 
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ratio P /P is known as a function of 0 for a given v . If 
mo o 

the atmosphere contains several components, then it is 
necessary to know P m /P Q f° r each molecular species as well 
as the proportion of each present. To take a simple example, 
if one considers a flat plate moving so that its surface is 
normal to the direction of motion, then 0=0 for the entire 
surface and we get 

P 

C D + 2(1 + J£) (Flat Plate). (2) 

o 

If the momentum of the reflected molecules is small compared 

to the incident momentum (P /P << 1) then C n - 2 , whereas 

if the molecules are reflected back along the direction of 

v q (specularly) with a speed equal to v q , then P^ = P q and 

Cp = 4. One, therefore, would expect the measured value of 

C D for a flat plate to be somewhere between these limits: 

2 < < 4. For a convex body it is conceivable that for a 

considerable fraction of the surface P is in the same 

m 

direction as P Q and hence is negative, leading to the 
possibility of values of less than 2. 

The present paper is concerned with the measurement in 
the laboratory of P /P as a function of v and 0 for 0 atoms 
incident on several surfaces, and the calculation of drag 
coefficients from the results of these measurements. These 
drag coefficients should be valid for the situation of a body 
moving through a rarefied, stationary gas of 0 for surface 
conditions equivalent to those in the experimental system. 
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SECTION II 
EXPERIMENTAL METHOD 


Most of the techniques employed in making the measurements 
have been described in previous reports, but they will be 
briefly given again here for completeness. 

The general procedure in measuring P m /P Q is to produce 
a beam of atoms or molecules having a known energy correspond- 
ing to satellite velocities (the energy in the case of 0 atoms 
is 4-10 eV), allow the particles to strike a test surface at 
a chosen angle of incidence, and measure the component of 
force on the test surface along the beam direction. If the 
rate at which the particles strike the surface (part. /sec.) is 
determined, then the force divided by the rate gives (P + P ), 
and since P q is already known from a knowledge of the energy 
and mass of the beam particles, then one has sufficient 

information to determine P /P . 

m o 


A. Equipment 

A schematic drawing of the entire experimental system 
is shown in Figure 1. The apparatus is mounted in two 
separate vacuum chambers, a test chamber and a beam chamber, 
which has the ion source, a focusing and deflection section 
and a mass analysis section attached to it. Mounted inside 
the beam chamber is another electrostatic focusing and 
deflection section along with either the neutralization cell 
and electrostatic collection plates or another focusing 
section, depending bn whether neutrals or ions were used to 
make a particular measurement. The test chamber contains the 
test surface which is mounted on a torsion balance used to 

test chamber is placed on a large concrete pier which is 
isolated from the laboratory floor to reduce mechanical 


4 



FIGURE 1. EXPERIMENTAL SYSTEM 



vibrations in the torsion balance. The two chambers are 
connected by a metal bellows which allows movement of the 
beam chamber so that the beam can be moved with respect to 
the test surface. The beam chamber is pumped continuously 
by a diffusion pump, but the principal pumping during operation 
is provided by a 6600 1/sec. cryopanel which is cooled to 
approximately 20°K by a Malaker Cryomite mechanical refrigerator, 

providing an operating pressure in the beam chamber of 0.8-4 x 

— 7 — 8 

10 Torr. Background pressure is approximately 2 x 10 Torr. 

The ion source is a magnetically confined oscillating 
electron bombardment type, patterned after the design of 
Carlston and Magnuson [6]. The pressure of the 0 ? gas in the 
ion source is in the low micron range and the anode current is 
typically around 10 ma at an anode voltage of 150 eV. The ion 
beam is extracted from the ion source through a 0.6-mm diameter 
hole by an extraction potential of 300 volts. 

After extraction the beam is focused and enters the mass 
analysis section, where the 0 + ions are separated from the 0 ^ + 
ions by a transverse magnetic field. A mass spectrograph of ^ 
ion collection current versus voltage across the electromagnet 
is shown in Figure 2. After more focusing and deflection in 
the beam chamber, the 0 + ion beam is decelerated to the proper 
energy and enters the neutralization cell which is located 
inside the bellows joining the two chambers. Within the 
neutralization cell, charge transfer takes place with the 
neutralizing gas krypton, and the beam leaving the cell is 
composed of 0 + ions that were not neutralized, neutralized 0 
atoms, and a small number of low-energy Kr atoms from the 
neutralizing gas. The charged particles are removed from the 
beam by a transverse electric field created by two parallel 
plates. The neutralizing efficiency of 0 + in Kr gas is only 
about 10% of that achieved with the N ? measurements and the 

X 

O' ion beam intensity is only about 15 to 20% of the ' beam 
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VOLTS ACROSS MASS ANALYZER MAGNET 
FIGURE 2. TYPICAL MASS SPECTROGRAPH 


intensity previously obtained. The resulting neutral 0 beam 
flux was such that reliable data could be gathered down to 
only about 2 0 to 25 eV energy. To achieve data down to the 
desired energies, the neutralization cell was replaced by 
focusing elements and the ion beam was utilized directly to 
perform the measurements. This will be explained in more 
detail later. 


The torsion balance which is used to measure the force 
produced by the beam on the test surface was patterned after 
one described by Pearson and Wadsworth [7] and uses electro- 
static damping and an optical lever [8] for measuring the 
angular deflection of the balance arm. This balance is 
relatively rugged but is capable of detecting forces as small 
as 2 x 10 dyne. The balance configuration used in our work 
is shown in Figure 3. The torsion fiber is 10-micron tungsten, 
and both damping plates are mounted on the same end of the 
balance arm. The test surface is mounted on the other end of 
the balance arm along with a momentum trap for measuring the 
beam flux. The procedure is to allow the beam to enter the 
momentum trap, measure the balance deflection, and calculate 
the beam flux under the assumption that the molecules leaving 
the trap have a Maxwellian velocity distribution which is 
characteristic of the temperature of the trap. (The results 
are not very sensitive to the precise validity of this assump- 
tion since the average momentum of nearly thermally accommodated 
molecules leaving the trap is much less than that of the incident 
molecules.) The beam is then moved upward mechanically so that 
it strikes the test surface and the balance deflection is again 
observed, giving the force produced by the beam. This method 
of measuring the beam flux has the advantage that the absolute 
calibration of the balance is not needed, since two balance 


average 


momentum transferred to the test surface by the molecules of 


the beam. In each measurement the ion beam is turned on and 
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off and the resulting balance deflection is taken as a measure 
of the force of interest. The effect of all forces on the 
balance other than that produced by the beam (such as that 
caused by the neutralizing gas) are thus eliminated. The 
momentum trap is 4.1-cm long with an apex angle of 22° and 
with a 0.9-cm diameter hole. The ratio of the hole area to the 
total internal surface area of the cone is 0.052, so that one 
would expect entering molecules to experience, on the average, 
around 20 collisions with the surface of the box before leaving. 
This means that if one assumes that the thermal accommodation 
coefficient for a single collision of a molecule with the inner 
surface of the cone is greater than around 0.2, the assumption 
that the molecules leaving the box have an average velocity 
that is characteristic of the box temperature is well satisfied, 
Although there have been no measurements of the accommodation 
coefficient in the eV energy range, the measurements at thermal 
energies [9] indicate that -for gas covered surfaces the 
accommodation coefficient is generally greater than 
Theoretical consideration of the particle-surface interaction 
as being hard-sphere at eV energies also leads one to think 
that the accommodation coefficient for these studies (with or 
without adsorbed gases) will be in excess of 0.5. All measure- 
ments described here were performed with test surfaces at room 

temperature . 

As stated earlier, the principal difficulty in using the 
electrical acceleration method, described in the previous^ 
section, is in achieving an adequate beam flux for detection 
by the torsion balance. One is concerned, then, with a signa - 
to-noise problem where the signal is the force produced by the 
beam on the test surface and the noise is that of the torsion 
balance. As an example of the magnitudes involved for a 10 e 
0 beam with ' a from , the 

ion current measured at the collection plates and an estimated 


10 


3% neutralization efficiency), the force produced on the test 

— 8 

surface would be about 7 x 10 dyne, under the assumption 
that the reflected molecules carry away negligible momentum. 

If the reflected molecules have appreciable momentum compared 
to the incident momentum, the force would then be larger. The 
minimum measured rms noise of this type of torsion balance is 
around 2 x 10 dyne [7], and this figure is just about equal 
to what one expects for the Brownian motion of the balance vane. 

The most important characteristic of the balance is its 
noise level expressed in force units. The period of the balance 
in its torsional mode is also a consideration since this deter- 
mines the length of time required for oscillations to be damped 
out after the beam is allowed to deflect the balance. It is 
also desirable to have the balance reasonably rugged so that a 
minimum amount of time is spent in construction and testing of 
each new balance. 

With a 10-micron tungsten torsion fiber the balance is 
sturdy enough to be constructed without special equipment and 

can be handled easily. A typical rms noise level in angular 

-7 

units for this balance is around 5 x 10 radians, or, expressed 

— 8 

in units of force, it is 4 x 10 dynes. At times of minimum 

external disturbance, the noise level of the balance was found 

— 8 

to be about equal to the Brownian limit of 2 x 10 dyne. The 
larger noise values result from pressure variations due to the 
irregular pumping characteristics of the diffusion pump and 
from vibrations reaching the balance. The torsional period of 
the balance depends on the moment of inertia of the vane 
assembly which is slightly different for each new balance but 
a typical value for the period is 15 seconds. 

In a measurement the beam is allowed to strike the test 
surface and the balance deflection is recorded. The beam 
chamber is then moved downward mechanically so that the beam 
enters the momentum trap, and the corresponding balance 
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deflection is recorded. The ratio of these two deflections 
then gives 

P + P 

R = ^- -p - , (3) 

o a 

where P is the momentum due to the atoms leaving the momentum 

3 

trap and has a maximum calculated value for these experiments 

of around 0.07 P . If one assumes that the molecules collide 

o 

with the walls of the momentum trap a sufficiently large number 
of times that they are in thermal equilibrium with the walls 
and leave with a corresponding velocity and angular distribu- 
tion, then by knowing the temperature of the walls, P can be 
calculated, and 


P 

p— can be easily computed from 
o 

P P 

= R(1 + —) - 1. (4) 

o o 


A straightforward investigation of the dependence of the 
ratio P m /P 0 on the angle of incidence would involve mounting 
the test surfaces on the balance vane so that the beam molecules 
strike the surface at the chosen angle , but with the beam 
direction still perpendicular to balance vane. This means that 
the balance must be modified or reconstructed for each new 
angle. At angles of incidence less than 30° this procedure was 
followed with success. At larger angles, however, it was found 
that the combination of the mass of the momentum trap and that 
of the larger test surface needed to intercept all of the; beam 
passing through the collimating aperture caused the balance to 
be intolerably noisy. Since all the surfaces studied^ gave \ 
essentially identical results it was decided that the measure- 
ments at larger angles would be performed by eliminating the 
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momentum trap and using one of the surfaces studied earlier 
(at 0 = 0°) as a reference surface, thereby reducing the mass 
mounted on the balance arm. The ratio P /P for the larger 
angles can thus be obtained by comparing deflections for 
the inclined surface and the reference surface, and then 
using tie results obtained previously which provided a com- 
parison of the reference surface with the momentum absorber. 
If the deflection for the inclined surface divided by the 
deflection for the reference surface is called S, then 


P + p 
o m 

P + P° 
o m 


1 + 


1 + 



(5) 


where P° is the value of P for the reference surface for 
m m 

0=0. Also, 


R 


o 


P + P° 
o m 

" P + P 
o a 



( 6 ) 


where R° is the value of R for the reference surface at 
0=0. Eliminating P°/P q between expressions for S and R° 
and solving for P m /P o one gets 


P P 

p— = S R° (1 + p“) - 1* (7) 

o o 

This is the expression normally used for obtaining P m /P Q for 
the larger angles of incidence where one is comparing the 
force on an -inclined surface to that on a reference surface. 
The above expressions were employed for all of the Nj measure- 
ments. 
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Regardless of whether atomic oxygen entering the momentum 
trap sticks, recombines or does not react with the surface, 
for an incident beam with energies over 4eV, P a /P Q is theoretically 
less than 0.07, and in this report P is assumed equal to 0. 

In the case that P /P is higher than zero, the corresponding drag 

3. O 

coefficient would not be increased by more than 3%. Hence, 
for a test surface mounted with the momentum trap, it is 
assumed that: 


P 

— = R - 1, (8) 

P 

o 

and for the larger angles of incidence where one is comparing 
the force on an inclined surface to that of a reference surface, 



S R° - 1. 


( 9 ) 


To perform the desired measurements in a reliable manner, 
the signal-to-noise ratio should be at least 10. As mentioned 
earlier, the force or the resulting deflection from a 10-eV 
0 beam is about three times the noise background of the torsion 
balance. Thus measurements were possible using the method 
utilized for the ^ results, and called the D.C. method for 
neutrals, down to about 75 eV. An A.C. or swinging method 
could be utilized down to about 20 eV before the balance noise 
and drift became too severe, and the D.C. method using 0 + ions 
rather than neutrals was employed down to 4 eV. These three 
techniques will now be described more thoroughly. 

-B-; B— C . M e-tfeod— 1 r ) 

The neutralization cell exit is located about 10 cm 

I 

from the test surface on the torsion balance. 
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This allows a considerable distance within which the 
neutral beam can diverge appreciably, especially at the low 
energies where the ion beam before neutralization is expected 
to be rather divergent. The diameter of the beam at the test 
surface at low energies was found to be around 1-1.5 cm. Since 
it is difficult to construct a satisfactory balance with the 
test surface and the entrance aperture of the momentum trap as 
large as this, some means was necessary to collimate the neutral 
beam before it reached the balance. This was accomplished by 
placing two 0.65-cm dia. collimating holes, one above the 
other, just before the balance. The diameter of the holes is 
such that all of the beam passing through the top hole will 
strike the test surface and all of the beam passing through 
the bottom hole will enter the momentum trap. The measurement 
is then performed by closing the bottom hole with a shutter, 
allowing a portion of the beam to pass through the upper hole 
and strike the test surface, and recording the corresponding 
balance deflection. The beam chamber is then moved downward 
a distance equal to the separation between the two holes, the 
top hole is closed and the bottom one opened, the same portion 
of the beam is allowed to pass through the bottom hole and 
enter the momentum trap, and the balance deflection is recorded. 
The ratio of these two deflections, then, gives the value of R. 

Because of the nature of the method used to produce the 
beam there are forces on the balance in addition to the desired 
force. This requires that the method used to obtain the balance 
deflection should eliminate any effects due to these extraneous 
forces. First, there is a force on the balance produced by Kr 
gas effusing from the neutralizing cell, which may be con- 
siderably larger than that caused by the particles of interest. 
Second, there can be a force produced by high-energy neutral 
particles that were produced by charge transfer of beam ions in 
the residual gas of the beam chamber at points within the 
electrostatic focusing system where the ion energy is higher 
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than the desired energy. One must then have a method of 
obtaining balance deflections which are due only to the 
desired neutral beam atoms and are not affected by the 
magnitude of these extraneous forces. This is accomplished 
by allowing all of these particles to strike the balance and 
then measuring the balance deflection that results when the 
particles of interest are prevented from reaching the balance. 
This in effect allows one to ignore the effect of the unwanted 
particles. The 0 atoms of interest (which in previous dis- 
cussion we have called the beam) are prevented from reaching 
the balance by changing the potential on an electrode just 
before the neutralization cell so that the ions cannot enter 
the cell. This eliminates the force on the balance caused by 
the neutral atoms formed by neutralization of these ions, but 
does not affect the forces due to the effusing gas molecules 
and the high-energy neutrals. The corresponding balance 
deflections caused by turning the ion beam on and off in 
this manner at 30-second intervals is then the desired 
deflection. This procedure is, of course, repeated for both 
the test surface and the momentum trap. 

Generally four or five measurements are taken with the 
beam striking the test surface , then a similar number with 
the beam entering the momentum trap, and then another set with 
the test surface. The average for the test surface is then 
compared to the average for the momentum trap. The fact that 
a complete measurement includes two sets for the test surface 
tends to minimize the effect of slowly changing beam conditions. 

C. Swinging Method (Neutrals) 

At beam energies of less than about 75 eV, the deflection 
of the tor sion balance was not sufficient when compared to the 
background noise deflections to make reliable mea'suremertts-j — so — 
an A.C. or swinging method was utilized. Briefly, the technique 
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used is to remove the damping from the torsion balance and 
turn the beam on and off in phase with the natural resonant 
frequency of the balance. In essence a constant force is 
applied in phase with a simple harmonic motion, thus causing 
the amplitude of the motion (rotation of the balance arm) to 
build up at a constant rate. A recorded trace of the amplifier 
output as this "in phase" cycling procedure is performed is 
shown in Figure 4 along with a typical trace of the D.C. method 
As before, a number of these "build-ups" are recorded with the 
beam striking the test surface, then a similar number with the 
beam entering the momentum trap , and then another set with the 
test surface. The slope of the plot of the magnitude of the 
rotation of the balance versus the number of swings then yields 

the average rate of the build-up. The slope is generally 

/ 

determined by a computer program using the least-squares method 
A ratio of the average of the slopes for the two beam positions 
then gives the desired value of R. 

Results using this technique compare within experimental 
error to results obtained using the D.C. method at the higher 
energies and reliable data is possible for neutral 0 atoms down 
to about 20 eV. 

D. D.C. Method (Ions) 

The mechanics or procedure for this method is identical 
to the D.C. method for neutrals, however advantage is taken of 
the fact that the ion beam intensity is much greater than the 
neutral beam intensity, thus permitting large signal-to-noise 
ratios down to less than 4 eV. The possibility of using ions 
rather than neutrals is based on the idea that positive ions 
are neutralized by an electron attracted from the surface 
immediately prior to impact on a conducting surface. The 
electron is attracted from the surface due to the electrostatic 
force when the ion is several angstroms away [10-11]. Our data 
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FIGURE 4. TYPICAL DATA TRACES 







taken at energies greater than 20 eV indicate excellent agree- 
ment with neutal data. However, a small energy correction of 
about 2 eV must be made to take into account the effect of the 
attraction of the ion due to the electrostatic image force prior 

to the neutralization. The ion beam, neutralized near the tar- 

3 6 

get surface, receives an additional energy, AE, where AE = — eV 
and d in Angstroms is the distance from the target. It is 
assumed that d is 1.8 8. 

When a measurement is made with ions instead of neutrals 
striking a test surface mounted at an angle, anouther phenomenon 
must be considered. As the ion approaches an angular suface , 
the electrostatic image force tends to bend the parth of the ion 
and cuase the trajectory to be more normal to the surface. The 
closer the ion comes, the greater the force and the more the 
trajectory changes, until the ion is neutralized several angstroms 
from the surface. The neutralized ion then continues on a straight 
line path to the surface. The amount of the trajectory change is, 
of course, dependent upon the angle of the surface and the energy 
of the incoming ion. As an example, a 4 eV 0 ion approaching 
a 75° surface, measured from the normal to the approach direction, 
is bent so that the actual impact occurs as though the surface 
was at 53°. From a measurement standpoint, this presents a 
difficulty, since the measured deflection of the torsion balance 
due to momentum exchange is not due to the desired component along 
the actual ion impact direction, but rather is that along the 
original ion beam direction. 

To circumvent this problem encountered with measurements 
using ions on angular surfaces, the following procedure is 
utilized. Data for curves of momentum transfer ratio versus 
energy at a given angle are taken for N 2 , N 2 + , and 0 + . For 
a given energy, the corrected impact angle, due to the bending 

a + + 

trajectory of the ion, of both 0 and N 2 are nearly the same, 
i.e., less than one-half a degree difference even for a 75° 
test surface. Assuming that the neutralized N 2 molecule and 
0 atom behave in a similar manner at the surface, i.e., the 
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reflected momentum patterns are similar for the same angle of 
incidence on the same surface for identical energies, then 
the ratio of measured S values for ^ to N£ + multiplied by 
the S value for 0 + yields the desired value of S for neutral 
0 for the desired angle of incidence of the test surface. This 
procedure is followed for all the non-zero angle of incidence 
ion results, and the validity of the above assumption is 
indicated by the excellent agreement of the ion and neutral 
data for the Echo I surface. 

To enhance the ion beam intensity, especially at the low 
energies, the neutralization cell was replaced by a lens and 
focusing section. This enabled the ions of desired energy to 
travel about 5 cm instead of 17 cm before striking the test 
surface and thus reduced the loss of beam intensity due to 
space charge spreading. Because of the changed angular 
geometry and to reduce the possibility of part of the beam 
missing Ihe target, the size of the collimating aperture at 
the shutter was lowered from 6.5-mm to 4-mm diameter. However, 
it was found that some of the beam did in fact miss the 
angular surface, so the torsion balance was electrically 
insulated to ena-ble the measurement of the ion current striking 
the surfaces. The balance deflections were then normalized to 
these currents to take into account the portion of the ion beam 
which did not strike the angular surface. 

E. Surface Condition 

The measurements were performed for several different test 
surfaces. The entire question of surface condition in experi- 
ments such as these involves a number of uncertainties. In 
considering the application of the measurements to satellite 
studies of the density of the earth's upper atmosphere, one 

would like to make the measurements~U's±ngr~safflpies— ©f— acixiaJ 

satellite surfaces which have the same surface condition as 
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that of the satellite in orbit (especially regarding adsorbed 
gases on the surface). There are two reasons why this 
desirability cannot be achieved at present. First, the 
condition of the satellite’s surface is to a large extent 
unknown. It depends on the preparation of the satellite, its 
environment in orbit, and possible continual emission of gases 
from portions of the satellite. Second, even the most advanced 
laboratory techniques are not presently capable of specifying 
precisely the condition of a surface under study. It is possible, 
however, that some aspects of the molecule-surface interaction are 
not especially sensitive to the exact nature of the surface, 
particularly aspects that involve averages over a number of 
parameters. Since the momentum transfer measurements described 
here provide a rather coarse study of the interaction, the 
following philosophy has been adopted with regard to surface 
condition. The measurements are performed for several test 
surfaces, but the exact condition of the surface is not rigidly 
controlled. The surfaces are handled carefully before placing 
them in the vacuum system so as to prevent their being con- 
taminated by oils, fingerprints, etc., but no attempt is made 
to remove adsorbed gases from the surfaces after they are in 
the vacuum system and the measurements are performed at 

_7 _ 0 

pressures (5 x 10 - 1 x 10 Torr) such that a clean surface 

(no adsorbed gases) cannot be maintained. If the results of 
the measurements indicate that only the gross character of the 
surface (such as surface roughness) affects the momentum transfer, 
then one might conclude that the surfaces can be adequately 
characterized for this particular type of measurement. Measure- 
ments that investigate finer details of the interaction, such as 
the angular and velocity distribution of the reflected particles, 
may require considerably more accurate surface characterization. 




\ 
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SECTION III 
RESULTS 


The objective in these measurements is to investigate 
P m /P G as a function of molecule energy and angle of incidence 
for two test surfaces. The surfaces used were samples of 
material used in the earth satellites Echo I and Explorer XXIV 
(no paint). Photomicrographs of these surfaces taken with a 
scanning electron microscope to show the gross roughness are 
shown in Figure 5. The results of the measurements are 
presented in Figures 6-8, in which the ratio P m /P Q is plotted 
versus the incident energy E . The measurements extend up to 
an energy of 200 eV since these results are rather easily 
obtained and they indicate the high energy limit of the 
momentum transfer. One of the principal factors in determining 
the nature of the particle-surface interaction is the ratio of 
the masses of the incident molecule and the surface atoms that 
it strikes. For this reason a measurement was previously made 
for a gold surface (mass number 197) at 0 = 0° to see if the 
results are affected by a large change in the mass number of 
the base material. The fact that the results for gold were 
essentially the same as for the other surfaces indicates that 
under the conditions of these measurements the interaction with 
adsorbed gases appears to predominate. 

In Figure 6 for the Echo I surface error bars are shown 
for the measurements at 0 = 0° and 0 = 75°, and are repre- 
sentative of the corresponding uncertainties in the other 
measurements. The error bars give the standard deviation in 

the mean value of P /P as calculated from a number of measure- 

m o 

ments at a given energy, and therefore represent the result of 
random fluctuations from the mea n value caused by system noise, 
etc. No inclusion has been made of possible systematic errors 
in the measurements, but it is felt that these should be small 
since the measurements involve a comparison of two determinations 
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FIGURE 6. MOMENTUM TRANSFER RESULTS FOR ECHO 





EXPLORER XXIV (NO PAINT! 
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FIGURE 7. MOMENTUM TRANSFER RESULTS FOR EXPLORER XXIV (LINES VERTICAL! 



ON EXPLORER XXIV (NO PAINT) 



FIGURE 8. MOMENTUM TRANSFER RESULTS FOR EXPLORER XXIV (LINES HORIZONTAL) 




of the same type of quantity (balance deflection) which were 
repeated a number of times for each surface and angle over a 
period of several months. 

It is seen in the photomicrograph of the Explorer XXIV 
surface in Figure 5 that there is a pattern of parallel grooves 
and ridges. It was found that for the measurements of P /P 
for this surface at the larger angles, the orientation of 
this pattern affected the results. The measurements shown in 
Figure 7 for 45° and 75° were taken with the grooves vertical 
(projection of the grooves on the normal to the beam direction 
is vertical), while those of Figure 8 were taken with the 
grooves horizontal. 

Drag Coefficients 

One of the purposes of these measurements has been to 
allow one to calculate drag coefficients for bodies moving in 
free-molecular flow with speeds in the satellite range. With 
the experimental results of the last section and Equation (1) 
one can calculate drag coefficients for a body of any convex 
shape. In this section such calculations for spherical bodies 
are performed. 

In order to do the integration indicated in Equation (1) 
over a sphere, it is necessary to know P m /P Q as a function of 
0. Figure 9 shows a typical plot of P m /P Q versus cos0 taken 
from Figure 6 for two energies. For a sphere. Equation (1) can 
be written as 

tt P 

0^ = 2[1 + 2 / T p^- sin0cos9d0]. (10) 

o o 

Taking the data points of Figure 9 and fitting them with 
straight-line segments (three or less) in the range 0.26 £ cos0 
<_ 1.00 one can easily calculate the contribution to the drag 
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coefficient for parts of the sphere where 0 is less than 75°. 

Since no data was taken for angles greater than 75° , then 

some extrapolation to larger angles must be used. Consider 

two such extrapolation procedures as bounds on the actual 

behavior of the curve in this region: (1) an extension of 

the straight-line segment used for angles slightly less than 

75° to angles 75° <_ 0 _< 90° will give the upper bound to the 

actual curve; and (2) a straight line drawn from the data 

point, at 75° to P /P = -1 at 0 = 90° will give the lower bound. 

mo s 

Plots of the drag coefficients obtained by these two methods are 
shown in Figures 10 and 11. It is seen that the two extremes 
in extrapolation amount to about a ±1 percent difference in the 
value of the drag coefficient. 

In Figure 12 the curves of drag coefficient for 0 versus 
energy have been plotted for the Echo I and Explorer XXIV 
surfaces, and for comparison, the ^ results on the Echo I 
surface are also shown. These represent averages of the two 
limiting extrapolation procedures, and for the case of 
Explorer XXIV, an average also of the two orientations of 
surface grooves (or lines). It is seen in Figure 12 that 
for a given energy (or velocity) the value of for the two 
surfaces vary by less than 2%. A large part of this variation 
is probably attributable to differences in gross surface 
roughness as observed in the microphotographs of Figure 5. 

When a surface like that of Explorer XXIV is mounted at some 
angle and has the lines or grooves oriented in a horizontal 
direction, then from a momentum exchange viewpoint, the angle 
effectively appears to be less than the measured angle. See 
Figure 13. The result would be a higher value of P m /P o when 
compared to a "smooth" surface or a surface with the lines 
oriented vertically. This effect was also noticed in the ^ 
momentum transfer results. 

The error bars in the low energy data (shown typically in 
Figure 6)would lead to about a ±2 percent uncertainty in the 
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FIGURE 10. DRAG COEFFICIENTS FOR ECHO 


ON EXPLORER XXIV (NO PAINT) 
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PARTICLE ENERGY 
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FIGURE 12. AVERAGE DRAG COEFFICIENTS 




FIGURE 13. SIDE VIEW OF SMOOTH SURFACE AND SURFACE WITH GROOVES 
HORIZONTAL 

drag coefficient, which along with the uncertainty in the large 
angle extrapolation procedure would lead to a total uncertainty 
in the drag coefficients of about ±3 percent for the lowest 
energies of the curves of Figure 12. 

Consider now the question of estimating drag coefficients 
for surfaces other than those for which measurements were 
described above. From the fact that the 0=0° results seem 
rather insensitive to both the base material of the surface 
and the gross surface roughness, the variation in value of 
drag coefficient seen for the surfaces in Figure 12 probably 
comes about mainly due to the different large angle behavior 
of the results. It is very likely that surface contour 
differences are what lead to this behavior, although the 
distance scale on which these contour differences are most 
important is not known. A good guess for an unknown surface 
might be arrived at by comparing a photomicrograph of the 
surface with those shown in this report and in the Nj report [5], 
and from the appearance of the surface roughness, estimate the 
drag coefficient by inspection of the measured results. 
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SECTION IV 
CONCLUSIONS 


The drag coefficients for the Echo I and Explorer XXIV 
spherical surfaces in an 0 environment have been experimentally 
determined over an energy range of 4 to 200 eV. The following 
results were noted: 

1. The drag coefficients studied are seen to be 
slightly greater than 2 at energies corresponding 
to satellite velocities. 

2. The value of the drag coefficients for 0 on the 
two surfaces vary by less than 2%, and much of 
this variation is probably attributable to effects 
of differences in gross surface roughness. 

" v t ' r . 

'3. The ' value of the drag coefficients for Nj and 0 on 
the Echo I surface differ by less than 1.5%. 

4. The value of = 2 for spheres obtained from both 
the and 0 experimental results comes about 
because the appreciable positive values of P m /P Q 
for small angles is to a large extent cancelled by 
the negative values at larger angles. A body of 

a different shape might therefore give values of 
Cq which are significantly different from 2. 

5. The validity of using ions instead of neutrals to 
increase the beam intensity for momentum measure- 
ments is justified by theory and by the excellent 
experimental results observed when both ions and 
neutrals were used on several different surfaces. 
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